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ABSTRACT 


An emalysls Is given of the errors in Hadamard 
spectroscopy that are caused by transparent silts in the 
mask being systematically wider or else narrower than 
they should be. It is shown that if the input spectrum 
consists of a single line, the distorted spectrum that is 
actually calculated consists of this line, plus four small 
blips. Wnen the transparent slits are too wide, these 
blj-ps are of equal height and the same sign, one pair 
surrounding the line, and another pair displaced a certain 
distance from it. When the slits are too narrow, the dis- 
placed blips have the same amplitude but are negative. 

Ihe response to an arbitrary input spectrum is then 
(istermlned from this. Ihe same method of analysis may 
also be used to handle other types of errors. 
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ERRORS IN HADAMARD SPECTROSCOPY OR IMAGING 
CAUSED BY IMPERFECT MASKS 


Introduction 

In the course of running a series of laboratory cali- 
bration spectra with a mercury vapor source at 1.7 microns, 
we observed a consistently occurring pair of negative dips, 
displaced a fixed number of spectral element positions to 
one side of the mercury vapor doublet Fig. 1(a). These 
looked like absorption features except that the intensity 
at these positions of our derived spectrum v'as actually 
negative rather than Just zero. 

Or. further euralysls, we found that moving the grating 
shifted both the spectral lines and these negative dips 
by equal amounts. However, when the spectral lines were 
moved to the extreme left, the two dips reappeared at the 
extreme right of the spectrum. This immediately suggested 
that we were dealing with a coding rather than with an 
optical error. The encoding mask has a cyclic code and a 
coding error could therefore be simply cyclic in this 
fashion, whereas stray light effects, for example, would not 
be as likely to produce a cycling error. 

A further piece of evidence, showing that light in- 
tensity did play an important role, was obtained when the 
mercury vapor doublet was moved so far to the left edge of 
the spectrum that one of the lines no longer was transmitted 
through the spectral mask. Wot only did that component of 
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the doublet then disappear from the spectrum, but the cor- 
responding negative dip also vanished. From this we concluded 
that the error was dependent both on the light intensity and 
OP our encoding scheme. 

This still left a number of possibilities open. 

After a number of attempts to find the correct explematlon, 
we finally thought of the following possibility. During 
tiie manufacture of the masks, whether by deposition of metal, 
or by removal of metal through an etching process, it is 
possible to obtain a systematic error that leaves each of 
the opaque portions of the mask either too wide or too 
narrow by a fixed amount. Essentially there is too much 
metal deposited, or too much metal etched away at each edge 
separating an open or a closed slit, and each of these edges 
is therefore displaced by a fixed amount. Independent of 
whether the open slit is wide or narrow. For example, if 
the open slits are too narrow, then the light passing through 
an open slit position may have Intensity values I^, I^(l-e) 
orIj^(l-2e) depending on whether this particular open position 
is bounded by two other open slits, has a closed slit 
bounding it on one side only, or bounds on closed slits on 
both sides. A computer simulation of this type of defect 
showed, in fact, that we had hit on the correct i.olution. 

Fig. 1(b) shows the computer response when provided with 
synthetic data expected for a single spectral line ajid a 
perfectly constructed 255-element S-matrlx encoded mask^. 

Fig. 1(c) shows the same output vdien the computer is pro- 


vlded with synthetic data corresponding to a mask in which 
each transmitting slot is too small by 0.1 slit widths 
at each edge. Fig, 1(d) shows the output for a similarly 
defective mask when the open slits are too wide by the 
same amount. 

While the source of the spectral errors had thus 

I 

been located, we felt the need for a more general treatment 
of this class of defect, since the displacement of the 
spectral dips presumably could be quite different for dif- 
fering encoding masks. We also wished to know what would 
happen if the transmitting slits were systematically too 
wide rather them too narrow. Finally, since image encoding 
masks cem make use of the same codes used in spectroscopy, 
we felt that a general mathematical analysis would also 
be germane to a wider range of optical problems outside of 
spectroscopy. 

The remaining sections of this paper therefore deal 
with this ema lysis. The results are remarkably simple. 
Independent of the particular S-matrix mask to be used there 
are always precisely four false blips in the response to 
any single element Impulse. The amplitude of these blips 
is always the same for a fixed narrowing or widening of the 
transmitting slits. Two of the blips always surround the 
main impulse and a pair of adjacent blips always is some 
distance removed from the impulse. We show how to compute 
the distance of this dirnlacement. The amplitude of the 
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displaced blips is positive when the transparent slits are 
too wide and is negative when the slits are too narrow. In 
contrast, the two blips surrounding the main impulse always 
are positive. 


(I ) Mathematical analysis if no errors occur in the mask 

We shall assume that there are n » - 1 mknowns. 

Let 

Xq, x^, . . . , denote the unknowns, 

^0*^1* * * **^n-l measurements, 

Zq, z^, . . ., denote the estimates of the unknowns. 

Also we write x « (Xq, . . . , 'i = \yQ> • • • » 

T 

z » (Zq, . . ., Zj^_2.) * where the T denotes transpose. These 
quantities are related as follows. First, is a Hadamard 
transform of x» given by 

Z = Sx (1) 

where S is a clrculant matrix whose general definition is 
given below in (III). As an example, if n = 15» a typical 
matrix S is shown in Fig. 2. The blank entries are zero. 
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Pig. 2 

An Example ol the Matrix S When n ■ 15 
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Equation (1) apsumcs that UiCrff lo no detector error, and 
that the mask is made correctly. Second, z is the Inverse 
Hadamard transfon.i of jr, given by 

z = S 


Of course in this case 

z_ = S”^Sx = X, 

and the unknowns have been determined exactly. 


r 
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(II) Errors caused when the transmitting silts are too wide 
Now suppose that the transparent slits in the 
mask (which correspond to the l*s In S) are slightly 
larger than they should be. Thus some radiation manages 
to get through slits which sliould be opaque, if they are 
adjacent to transparent slits. More precisely, consider 
an arbitrary row of S, say 

... 11001011000011 ... . 

We assume that 

every 0 which is adjacent to two O's is unchanged, 

every 0 which is adjacent to a single 1 is replaced by c » & 

every 0 which is adjacent to two l*s is replaced by 2e . 

Here € is a small positive number which measures the fraction 
of radiation which manages to get through a slit which should 
be opaque and which is adjacent to exactly one transparent 
slit. Thus the preceding row would be changed to 

... 11c c 12c 11c 00c 11 ••• • 

Let S* denote the distorted matrix obtained from S in this 
way. Figure shows the distorted matrix corresponding 
to Fig. 2.. 
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Fig. 3 

The Distorted Matrix S* When n » 13 
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Now wo actually measure (instead of Eq. (1)) 

V “ S»x (2) 

and the estimate of the unknowns is 

z* = = S'^S*x. 

We wish to discover how z* differs from x. 

In fact it is enough to detennlnc z* when x has 
n single nonzero component. For suppose tlio input 

X = x^^^ = (1,0,0, . . .,0)'^ 


produces the response 



(say). 


(3) 


- t - 

z* - 

We call the Impulse response . Then since e'^ery row of 

S (and S* ) is a cyclic shift to the left of the previous 
row, it follows that 

X ■ (0,1,0, •t.,0) 

produces the response 

z* a 

and 

X « (0,0,1, . . .,0)'^ 

produces 

and so on. Here £ denotes a cyclic shift of one place 
to the left. Then an arbitrary input 

produces the response 

z* a + ... ♦ . (5) 

M ' 

Therefore, once t’le impuJse response z' ^ is Known, Uie 
response to an arbitrary input (4) is determined (by (5)). 
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Befcre giving the analysis of the 
let us calculate for the case n ■ 15. 

S"^, where - stands for - 1. 

Fig. 4 

The Matrix S~^ When n ■ 15 
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Then Is found from Eq. (5) to bo 

» (l-t,^,0,|€,^,0,0,0,0,0,0,0,0,0,^ ). (6) 

(III) The matrix S 

The theory of the matrix S Is Intimately connected 
with that of finite fields. However, we shall give a simple, 
self-contained description in several steps which assumes 
no previous knowledge of finite fields. The definition of S 
given here Is In fact the same as that given In references 
1, 2 and 3. 
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general case. 
Figure shows 
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ci) Choose a primitive Irreduciole polynomial . To begin 
with, we need to choose a primitive Irreducible 
polynomial 

p(x) ■ 1 + a^^x -f SgX^ -f . , , ■♦■X*”, a^ ■ 0 or 1, 

of degree m, where n - 2®-l, The definition of such a 
polynomial need not be stated here. For convenience 
Fig, 5 gives a short table of such polynomials 


Fig. 5 

Primitive Irreducible Polynomials 
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Extensive tables of such polynomials are available in the 
U 6 7 8 

literature * ■** * * , For our example with n ■ 15, we 
take p(x) a l«xix^. 

(11) Form a shift register . Next, w»> form a linear feedback 
shift register with m otages, whose? feedback cc^nnections are 
described by p(x). For example, the shift register 
corresponding to p(x) * lixtx^ is shown In Fig, 6, Here -p]- 
is a delay element which holds a 0 or a 1, and ® is a mod 2 
adder, ^ 
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Fig. 6 

Shift Rggloter Correoponding to 



(111) Fom the matrix M . We loa^ this shift register 
with the Initial state (1, 0, 0, . . . , 0) and let it cycle. 
Becaase p(x) wan chosen to he a primitive irn duclble 
polynomial, the shift register will go through 2”'-l states 
before repeating. It has period 2™-l. During one period, 
the successive states of the register give all 2™-l nonzero 
m-tuples. 


For example, the first few states of the shift 
register of Fig. 6 are shown In Fig. 7 . 

Fig. 7 

Successive States of the Shift Register of Fig. 6 
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We form an m x (2'"-l) matrix M whoae columns 
are the successive states of this shift register. For 
n - 13, for example, we find from Fig. 6*7 that 

0 1 2 3 4 5 6 7 8 9 10 11 12 14 


M - 
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0 0 0 ? 0 

0 0 1)0 
0 10 0 1 

1 0 0 0 1 
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In general we have 
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The columns of M consist of all 2 -1 distinct nonzero binary 

m-tuples. The rows are labeled shown. 

2 

Note that in Eq. (7)p Vg ■ * C and 

V4 - 


where 0 denotes componentwise addition modulo 2, without 
carries. Also v^^ = » f."^Vj^, where R denotes a cyclic 

shift of one place to the right. Therefore 


- 


o ^2 O « 0, 

or 

(r^O^ 01)^1 - 0. (9) 

In the general case, the same argum'nit applied to Eq, (8) 
shows that 


© • • • O O V >''1 - (1°) 

l.e. , 


P(r)v2 - 0 (11) 

(iv) The connection with Galois fields . (This section 

car be omitted by readers only Interested In the statement 

of the result and not In Its derivation.) The successive 

contents of the shift register in fact give all nonzero 

elements of the Galois field GF(2™). These elements can 

be written in two ways, either ar successive powers 

2 2"'-2 
l»a,a » . • • «a 

of a, where a is a zero of the polynomial p(x), and 
2*"-l 

Tt ■ = Ij or equivalently as a polynomial in a of degree 

less than m, by reading off the contents of the shift 
register. In the example with m «= 4 we have 
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Thus M Is actually a logarithm table for the field 
OF (2*")! E.g, column 7 of Eq. (7) says that 

1 + a + f 


or 

log^(l + a + a^) =* 7. 

(v) Definition of the matrix S . We can now define S, 
which is the (2”'-l) x (2"’-l) matrix whose first row is 
(the first row of M), and each successive row is a 
cyclic left shift of the previous row by one place. If 
denotes the i " row of S (calling the top row the 0 
row), we have 




for i = 0, 


2"’ -2. 


Indeed Pig. 2 is obtained In this way from Eq. (7). 

It con be shown^ that 

^ (2S-J), (12) 

where J is the all-ones matrix. In fact by adding a top_ 
row and a left column of O's to S and then changing I's to 

-I's and O's to +l's, we obtain a Hadamard matrix H 

2 -1 
satisfying H » (n+l)I. Figure 4 shows S when n ■ 15. 

(vl) Labeling the rows of S . It turns ouV^®, and this 

is a key point, that the rows of S consist of all linear 

combinations of the rows of M, with coefficients 

0 or 1 but not all 0. I.e., the rows of S are equal to the 

vectors 

Vi» . . 0 0 . . .,Vj 0 Vg O • • • 0% 

arranged in some order. The rows of S can therefore be 
labeled with the corresponding sum of the v^^'s. For example, 
the rows of S are labeled in this way in Fig. 2. We also 
use the same labeling for the components of £. Thus if 
the i^^ row of S is labeled 0 0 . . . 0v^, the i^^^ 

component of z is labeled in the came way. 

(IV) The impulse response 

We can now specify the Impulse response. If 
x^^^ a (1,0, ...,0)^ then the Impulse response 


is given by 


16 - 



( 1 “€ » Of » • » fOf'^f ■^f Of ,, m 



( 13 ) 


where the second and third components ^ occur at the 
coordinates which are labeled vj + and + Vg. In 
words, has the dominant component 1-c (where had 

the component 1), surrounded by two blips of height plus 
a second pair of blips of the same height which are displaced 
a certain distance. The proof of this will be given in section 
V. See Fig. 8, which shows z when n » 15 (in agreement with 


Eq. (6)). 


Fig. 8 

Impulse response when n = 15 
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To discover exactly how far the second pair of 
blips is displaced, we may proceed as follows. Because 
of Eq. (10) & (11), there is an isomorphism^^ between the 
rows of S and the nonzero elements of the field GF(2*”), 
as follows ; 

row Corresponding el<^ment of GF(2"*) 

Pq » 1 

pj^ « Vg » <e-^ a 

<= for 0 ^ i < 2*"-2 

Pm “ ^ 

2”'-2 ^ 

Therefore the second and third blips in the impulse response 
occur at the i^^ and (l4l)®^ coordinates of (starting 

the count at 0), where 1 is the solution of 

1 ©a"^ = (14) 

in the field 0F(2"'). Indeed, for this value of 1, we 
have 


■'l ©''m ■ ''l©^‘S “ ■‘‘S ' Pi- 

Fortunately the columns of M provide the Information needed 
to solve (l4). For example, when n = 15, we have 




- lb 


so the solution is i ■ 3« Therefore the blips occur in 
the 3***^ and coordinates of (calling the left-most 

coordinate the 0^^), in agreement with Eq. (6) and Fig. 8. 
The case n ■ 255 

S matrices of size 255x255 have been used in 
Ref. 2 and 3» corresponding to the polynomial 

p(x) » 1+x^+x^+x^+x®. 


The matrix M obtained from this is an 3x255 matrix, and is 
shown in the appendix (turned on its side and cut into six 
pieces). Using this appendix Eq (l4) becomes 



Therefore the irapulse response is 

= (l-e,|,228 0's,|,|,22 0's,|) 

(V) D erivation of the impulse respons e 

If x^^^ = (1,0, ...,0)^ then from Eq. (2), since 
S* is symmetric, the transpose of is 

= + €w, (say). 


where w Is a vector of O's, I's and 2's, and + denotes real 
vector addition. Observe that w Is given by 

w - Vj^ » Tvj + Vj^ * Rvj^, 

whe»*^ the bar denotes the complement of a vector and * stands 
for componentwise multiplication, 

' ^1 * Vg + . v^ 

» Wi + Wg (say), 

where w^ « v^^ * Vg* Wg ■ v^^ * v^^^. Then 

« S"^X* = S"^(v]^+€wJ’-few2) 

« S“^v^ i- e2^""'(2S-J) (w^4-wg) 

« + €2^”"'(GwJ+GWg) 

where G = 2S-J is a matrix of +l's and -I's, If p Is a 
typical row of S, let 

=<PiWj^>= pw^, Xg =<p,Wg>= pWg 

12 

denote the usual real inner products. Then the component 
T T 

of GWj^ + GWg corresponding to p is 


from (12) 

(15) 


20 


To show this, let a be an arbitrary row of 0, and let 
£. “ i) corresponding row of S, where jL denotes 

a vector of all I's. Also let w be an arbitrary vector 
of O's and I's, containing say h I's. Suppose there are 
1 coordinates where £ and w are both 0, J coordinates where 
£ is 1 emd w Is 0, and k coordinates where £ Is 0 and w 
is 1, as shown in the following picture. 

» a 

“ £ 

“ w 

T 

Hien aw « 1-k. But 1 = 1+k » h, so 

aw^ « 2<^,^ - h . (16) 

T T 

This Implies that the component of GWj^ 4 Gw^ corresponding 
to £ is 

2(^1 + ^ 2 ) ‘ 2'"’^ • 

The values of the inner products X, Xg and, from 
Eq* (15)» of the components of are easily calculated, 

and are shown in Fig. 9 * 
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nils completes the proof that the impulse response 
Is given by Eq. (15). 









VI. Errors caused by the transparent silts being too narrow 


Now let us consider the case where the transparent 
slits are too narrow. More precisely, we assume that in 
each row of S, 

every 1 which is adjacent to two I's is unchanged, 
every 1 which is adjacent to a single 1 is replaced 
by 1 - € , and 

every 1 which is adjacent to two O's is replaced by 

1 - 2 €, 

# « 

while the O's are unchanged. Let S be the distorted 
matrix’ obtained from S in this way. Fig. 10 shows the 
distorted matrix corresponding to Fig. 2, where a denotes 
1 - € 8 Uid b denotes 1 - 2 c. 

Figure 10 

The Distorted Matrix S Wlien n t.- 15 
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Now the Input 

X - - (1,0,0,. ..,0)'^ 

produces as output the Impulse response 
■ S S X ' ' . 

This Is given by 

■ (l-€, 0, ..., 0, - --J, 0, 0, (18) 

where the components ^ occur at the coordinates which are 
labeled euid + Vg. Fig. 11 shows when n ■ 15, 

Note that is the same as except for the sign of 

the second pair of blips. 

Figure 11 
( 2 ) 

Impulse response z' ' when n » 15 



e/2 



TTie proof of Eq. (18) is similar to the proof of 
Eq. except that instead of Eq. (13) one has 

2 ( 2 ) , ^(1) _ (2GVj^^ - GWj'^ - Gw^^) , 

where * Vg and * v^» Also Instead of 

Eq, (17)i one has that the coMponent of 2GVj^^ - Gw^^ - Gw 
corresponding to £ is 

^<£,vp - 2 <£,w^ - 2<£,w^> ~ 2'"'^ . 

The details of the-proof are left to the reader. 

/ 

The part of this study tnat was done at Cornell 
University was supported by WASA contract NGR 33-010-210. 



REFKRLNCES 


1. N. J. A. Sloanc, T. Fine, P. G. Phillips, and 

M. Harwlt, Appl. Optics, 8 , 2103 (I 969 ). 

2. J. A. Decker, Jr., Appl. Opilcs, 1^ 510 ( 1971 ). 

3 . M. H. Tal, D. A. Brlotta, Jr., N. Knmath, and M. Harwlt, 

A practical multi -spectrum Hadamard transform 
spectrometer, 

U. W. W. Peterson and E. J. Weldcn, Jr., Error-Correct inn 

Codes , 2. id Edition, (M.I.T. Press, Cambridge, Mass. 1972). 

5 . J. D. Alanen end D. E. Knuth, Sankya, Series A, 26, 

305 (1964). “7 

6 . R. W. Marsh, Tables of Irreducible Polynomials over 
GF(2) Through Degree I 9 (Office of Technical Services, 

Dept, of Commerce, Wash., D.C., October -24, 1957). 

7 . S. Mosslge, Math. Comp., 26, IOO 7 ( 1972 ). 

3. W. Stahnke, Math. Comp., 2J, 977 (1973). 

9 . Strictly speaking this Is a half-odder, since there 

are no carries. For more about shift registers, see for 
example Ch. 7 of Ref. 4. 

10. This Is a standard result In coding theory, and expresses 
the fact that the rows v,,...,v of M generate a 

cyclic co<le C, of lengthen dimension m, and 

Hamming distance 2*’'"! between codewords. The rows of 
S ore exactly the nonzero codewords of C. Gee for example 
Ref. 4, or the forthcoming book F. J. MacWllliams and 

N. J. A. Sloane, Combinatorial Coding Theory (North-Holland, 
Amsterdam, in preparation). 

11. In other words the code C Is isomorphic to 0F(2”'). 

] 2 . In fact we are here finding the weight distribution 
of the cosets of C which contain w, and Wp. For more 
about the v/elght distribution of tnese cosets see 
N. J. A. Sloane and R. J. Dick, IEEE Internat. Conf. 
on Commun., Montreal 1971# Vol. . pp. 36-2 to 
36-6 (IEEE Press, N.Y., 1972)? and E. R. Berieknmp and 
L. R. Welcli, IEEE Trans. Info. Theory, IT-lcl , 203 
(lv72). Another application of these cosets In studying 
Hadamard transforms Is given by E. R. Bcrlekamp, Bell 
Syst. Tech. J., 4^* 969 (1970). 


Figure 1 Caption 

Spectrum of the 1.7 micron mercury vapor doublet (a). Peaks 
extending downward from the top of the picture are negative 
values which ordinarily would fall the same amount below the 
abscissa, but, In our form of display, appear at the top of 
the graph, (b) Shows the response we would obtain to a single 
spectral line with a perfect mask, (c) Shows the response for 
a single line with the radiation simulated as passing through 
a mask with slits too narrow because each opaque mask element 
protrudes into the adjacent transparent slot by a tenth of a 
slot width. (d) Shows the effect of simulating silts that 
are systematically too wide. Note that the main spectral line 
has been placed In different positions for the synthetic runs 
(b), (c) and (d). 
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Al’I’KNDIX: llie (Jaloln Field GF(2®), defined by i- 0 
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•nie Galois Field GF(2 ), continued 










